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1 . Introduction 

Recently, spatial deformations of interaction couplings in low-dimensional quantum 
many-body systems have attracted much attention [T],[2l[3llll[SlEl[7llSlini [IHl [HI 
A smooth cutoff of the coupling encoded in the real space is generally 
expected to suppress the scattering at the boundaries [12]. Among the various schemes, 
the sine-square deformation (SSD) in one-dimensional (ID) quantum critical systems 
is of particular interest, because the many-body ground state of a system with SSD 
almost coincides with that of the uniform system with periodic boundary conditions 
(PBC). This is due to the nearly exact cancellation of the boundary scattering in the 
system with SSD nna. In a class of ID critical systems which are reducible to free 
fermions, moreover, the ground-state correspondence can be made exact. In fact, it was 
proved that the ground state of the spin-1/2 XY chain with SSD is identical to that of 
the uniform chain with PBC [1]. The same holds true for the tight-binding chains [12] 
and transverse held Ising chains at criticality [5]. More generally, the ground-state 
correspondence between the SSD and uniform systems can be explained in terms of 
conformal held theory (CFT), where the Virasoro generator L±i identihed with the chiral 
part of the Hamiltonian of the SSD system im annihilates the CFT vacuum [5]. Another 
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interesting aspect of the SSD is that it has a variety of applications; For instance, the 
nature of the quasi-localized edge excitations originating from the SSD is essential in a 
‘grand-canonical’ approach for hnite-size systems, where smooth magnetization curves 
of quantum spin systems can be obtained with very high accuracy Besides, very 

recently, a further relation between the SSD and CFT, and its relevance to string theory 
is also pointed out in Refs. mm 

A key to understanding the ground-state correspondence between the SSD and 
uniform systems is that each chiral part of the SSD Hamiltonian annihilates the many- 
body ground state of the uniform system with PBC [12]. Here, we remark that the 
plane-wave representation of the single-particle basis, which is a suitable description 
for the uniform system with PBC, plays a very essential role in verifying the above 
statement. However, as shown numerically in Ref. piaiH], low-lying single-particle 
states in the SSD system seem to be localized around the center or edges of the chain, 
which would be rather puzzling. A natural question is: “What is the unihed view of 
the physics behind the plane-wave description and the single-particle eigenstates of the 
SSD?” In this paper, we explore this question by examining the SSD of the free-fermion 
system dehned in ID continuous space. To this end, we formulate the SSD problem 
as the inverse problem of the quantum mechanical problem of a particle in the inverse 
sine-square (1/sin^) potential. An elegant way to solve the direct problem is to use 
supersymmetric (SUSY) quantum mechanics with the shape invariance [IHlIin]. With 
this, we show that a set of the single-particle eigenstates of the SSD problem can be 
obtained in closed form for specihc values of the parameter /i, which will be introduced 
in the next section. Then we show that the Slater determinant obtained by hlling 
these single-particle states up to the zero energy is identical to the Slater determinant 
corresponding to the ground state of the uniform system with PBC. Moreover, we discuss 
the relation of the SSD in a lattice with that in the continuous space, and clarify the 
role of the geometric correction involved in the single-particle wavefunctions of the 
SSD. Then, we demonstrate that a number of properties of the lattice SSD found in the 
previous study [8] can be well explained on the basis of the SSD wavefunctions obtained 
with the use of SUSY quantum mechanics. 

The rest of the paper is organized as follows. In section 2, we dehne the SSD 
problem of free fermions in ID continuous space. In section 3, we formulate the single¬ 
particle problem with SSD as the inverse problem of the 1/ sin^ potential problem. Then 
we obtain the single-particle eigenstates for specihc values of the chemical potentials. 
In section 4, we prove the ground-state correspondence between the uniform and SSD 
systems by showing that the single-particle eigenstates for the SSD can be expressed 
as linear combinations of plane waves. In section 5, we briehy comment on the single¬ 
particle states with positive energy in the SSD system. In section 6, we compare the 
SSD in the continuous space with that in a lattice, and discuss the role of the geometric 
correction. In the last section, we summarize the results obtained and discuss the 
prospects for future researches. 
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2. Definition of the SSD problem in a one-dimensional ring 


We start with the definition of the SSD problem in ID continnous space. Consider 
N fermions on a ring of length L. We denote by c^{x) and c{x) the fermionic field 
operators at position x {0 < x < L). They satisfy the canonical anticommntation 
relations: {c(a:), c'*'(|/)} = 6{x — y), {c{x),c{y)} = {c\x),c\y)} = 0. In analogy with 
lattice systems, the Hamiltonian for the system with a generic deformation is defined as 




dx c^{x) f{x) 


dx'^ 


/i c{x) 


( 1 ) 


where y is the “chemical potential” introdnced for later convenience, and f{x) encodes 
the spatially varying conplings. Here we set h = 2m = 1. If f{x) = 1, the systems is 
nothing bnt the nniform ring with non-relativistic free fermions. For the case of SSD, 
f{x) is defined as 

f{x) := sin^ (^^x^ . (2) 

Thus, the energy scale of the system smoothly approaches zero as x —)■ 0 and x ^ L, 
implying that the ring is effectively disconnected at a; = 0 (mod L). In this sense, the 
boundary conditions of the SSD problem can be interpreted as open ones. 

A general Wparticle eigenstate of T^cont is of the form: 


1$) = / dxi--- dxNc\xi)---c\xN)\0)^{xi,...,XN), 
Jo Jo 


( 3 ) 


where |0) is the state with no fermions. Since the Hamiltonian is quadratic in c 
and ci, the Wbody wavefunction in the coordinate representation $(xi, ...,xn) can be 
constructed as a determinant of single-particle wavefunctions (p{x) which are determined 
by the Schrodinger-like equation: 


- fix) ( ^ + /^ ) ^(x) = i^ix), 


( 4 ) 


where ^ denotes the eigenvalue. Note that the Hamiltonian here is not self-adjoint in the 
standard sense and the appropriate inner product for the SSD problem will be introduced 
in section 5. For later convenience, we assume either periodic or antiperiodic boundary 
condition for 9 ?(x), depending on the total number of fermions: ip{x + L) = (—l)'^<y5(a;). 
We denote by and the sets of single-particle eigenfunctions of the 

uniform (/(x) = 1) and SSD (/(x) = sin^ ^x) systems, respectively. For the uniform 
system, {v5jj“(x)} is of course described by the plane waves and is easy to obtain. For 
periodic boundary conditions, one finds that the linear space of is spanned by 


2mi 
exp ±i——X 


and for antiperiodic boundary conditions, 

( X2n + l)Ti 

exp ±i ---X 


(n = 0,1,2,...), 


(n = 0,1,2,...). 


( 5 ) 


( 6 ) 
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In contrast, to obtain for the SSD problem is an intrigning problem, because 

the translational symmetry is explicitly broken by the position-dependent factor f{x). 
However, for specihc values of the parameter /i, one can obtain them in closed form by 
exploiting the hidden SUSY in the problem. 


3. Single-particle eigenstates of the SSD problem 


In this section, we present a systematic procedure to construct single-particle eigenstates 
of the SSD problem. To this end, we introduce the following parametrization of 

?=-0UW-l). (7) 

and rewrite Eq. (jl]) with f{x) = sin^ jx as 

H^ifix) = nip{x), ( 8 ) 


where the “Hamiltonian” is dehned as 


(P /7r\2 /j(/3-l) 

^ dx"^ \L) sin^ 


(9) 


In the standard setup, “energy eigenvalue” fi is determined for a given parameter ^ (or, 
equivalently, /5), where the corresponding eigenfunction has a wavepacket-like shape. 
On the other hand, however, what we need here is to determine the energy eigenvalue ^ 
for a given chemical potential /r. In this sense, the SSD problem can be formulated as 
the inverse problem of the 1/sin^ potential problem. 

It is known that the 1 / sin^ potential problem can be solved algebraically since 
the system has SUSY and shape invariance [IH]. This property is also deeply related 
to the integrable many-body system with an inverse sine-square interaction solved by 
Sutherland [201 EH- Here, we briefly summarize the SUSY quantum mechanics for the 
1 / sin^ potential problem. The eigenvalue problem of Eq. ([8]) can be solved algebraically 
for non-negative integers (3 = 0,1,2,3---. Let us start with dehning intertwiners 

with the superpotential Wjs^x) := —/3 In (sin . Using these operators, the 

Hamiltonian in Eq. ([8]) can be written as 



According to SUSY quantum mechanics, it is well-known that and form the 

SUSY-partner Hamiltonians, and the shape invariance ensures the recursive relation 

+ (—( \{fi+ 1 )^ - fi'^] — H^+i - 


( 12 ) 
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In addition, it is easily confirmed that A^'ip^{x) = 0 with := (sin^x)^, which 

corresponds to the eigenfnnction of with the lowest eigenvalue. Note that Eq. OUl 
implies that the eigenvalue of is bounded from below by ft follows from Eq. 

(na that, for a given non-negative integer /3, the f'th eigenstate of is constructed as 

The explicit form of the eigenfunction can be written as 

cf (cos 

with the coefficient 

, ^ A. (2/3 + 2£-l)!(/j-l)! 

^ \ l) 2^(2/3 + £-1)!(/3+ £-!)!’ 

where {z) denotes the Gegenbauer polynomial of the £th order [221 [21]. The coefficient 
A/? can be obtained by demanding that A^'i/j^^{x) = (^) £{£ + 2(3)^jJ^^_^\{x) and the 
following recursion relation for the Gegenbauer polynomials are consistent: 

Tcf(a)=2/3CS(a). (16) 

Then, the corresponding eigenvalue of is obtained as 

h= (0'(/5 + ^)', (17) 

for £ = 0,1, 2, • • •. Thus, if the sum of (3 and £ is constant, say (3 -\- £ = q, then fi takes 
the same value among '0°(x), • • • , 'ipp~^{x), ■ ■ ■ ,'iPq{x). 

Here, we list in Table I the first few eigenfunctions of Hp with /3 = 0,1,..., 5, from 
which the relation between the 1/sin^-potential problem and the SSD problem can be 
seen clearly. On the one hand, one finds the eigenfunctions of Hj^ in the column labeled 
by (3. On the other hand, the functions in the row labeled by q give the solutions of the 
inverse problem (up to unimportant overall factors). For a particular chemical potential 
parameterized by /i = (f)^q'^, we obtain them as {£3 = 0,1,..., g) with the 

corresponding energy eigenvalues ^ = — (f)^/d(/d — !)• Note that ^g(x) can be thought 
of as the single-particle ground state of the SSD problem, while '^o(^) ^cid 'tp'£~^(x) are 
the plane-wave solutions with zero energy. 


(13) 

(14) 

(15) 




Table 1. List of eigenfunctions of Hp up to /3 = 5. P and AP in the B.C. column, respectively, indicate the periodic and 
antiperiodic boundary condition imposed on the eigenfunctions. The columns for /3 = 0 and 1 correspond to the plane-wave solutions 
of the uniform system. The functions in each row give the eigenfunctions of the SSD problem. For fixed g, they are identical to 
• • • , (x), • • • , fpQix) up to unimportant overall factors. 


02 
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In order to precisely see the correspondence between the solutions of the 1 / sin^- 
potential and the SSD problems, we further comment on the boundary conditions. In 
the literature, it is often assumed that /3 > 1 in Eq. (IH]). Here, we have added the SUSY 
partners for /3 = 0 and 1. For both /3 = 0 and 1, there is no potential term in Hp and 
thus the system is uniform in 0 < a: < L. However, an important difference between 
them is that the boundary conditions for /3 = 0 are Neumann, i.e., <yC^(0) = ^'{L) = 0, 
in contrast to the Dirichret ones for (3 = 1. Thus, the eigenfunctions for /? = 0 and 1 
are described by cosine and sine functions, respectively. In terms of SUSY, the lowest 
eigenstate of Hq is a singlet with zero energy, because it is annihilated by both Aq and 
Hq. Any other eigenstate of Hq has the corresponding eigenstate of i7i, which give rise to 
a doublet of eigenstates with the same eigenvalue. Thus, we can express the plane-wave 
solutions for a ring of length L by superposing the eigenfunctions of Hq and Hi. We 
recall here that we have imposed periodic (antiperiodic) boundary conditions on (p(x) 
depending on the total number of fermions N. Therefore, in Table I we should take the 
eigenfunctions in the row labeled by even (odd) q when N is even (odd). We also note 
that the boundary conditions for the eigenfunctions are indicated by P (periodic) and 
AP (antiperiodic) in Table I. In contrast to the cases of /3 = 0 and 1, the eigenfunctions 
flTT)) with (3 > 2 enjoy (p(0) = (^(T) = V^^O) = <^'(T) = 0, which reflects the fact that 
they have wavepacket-like shapes and are disconnected at x = 0 (mod L). 


4. Equivalence of two Slater determinants 


In this section, we give a proof that the Slater determinant of the plane waves, 
corresponding to the ground state of the uniform system with PBC, is identical to the 
Slater determinant of which are the eigenfunctions of the SSD problem. This 

can be thought of as a continuous version of the ground-state correspondence proved 
in a class of ID lattice models reducible to free fermions mum- In the following, 
we focus on the case where the total fermion number N is even. It is, however, very 
straightforward to extend the proof to the case of odd N. 

Let us hrst consider the many-body ground state of N (even) fermions on the 
uniform periodic ring with /i = (^)^Y^, which is the Slater determinant obtained by 
hlling the plane-wave states up to .^ = 0. It is easy to see from Table I that the plane- 
wave states with even q can be obtained as linear combinations of eigenfunctions in the 
columns labeled by /? = 0 and 1, which satisfy the PBC. Since fi = (^)^Y^, one can 
easily conhrm that the states in the q = N row with = 0 and 1 are just located at 
.^ = 0. Thus, the many-body ground state is two-fold degenerate, because either one of 
the two ^ = 0 states is occupied by a fermion. In the following, we assume that /? = 1 
state is occupied. Then, the single-particle states involved in the Slater determinant is 
written as 


^7\x) = 


cos 


sm 


(2n—A^H-2)7r 


(u = 0,l,...,f-1), 


(n = f,f+ 1,...,Y-1), 


X 


(18) 



where cos does not appear because we have assumed that /3 = 0 state is empty. 
Note that, if we consider the case of odd iV, the boundary conditions are antiperiodic, 
where we should take the wavefunctions only from q =odd sectors. In terms of 
the iV-fermion ground-state wavefunction is written as 

...,xn) = det • (19) 

1<2 J<iV 

We next consider a many-body ground state of the SSD problem with the same 
setup, i.e., N is even and n = As shown in the previous section, the single¬ 

particle eigenfunctions of the SSD problem is given by 

^f{x) = (/3 = 0,1, 2,..., iV), (20) 

where 'ijj^~‘^{x) is dehned in Eq. ([13]). The many-body ground state is then described 
by the Slater determinant of (1 < /9 < -^)) which span the subspace of 

Here, we again assume that not /? = 0, but /9 = 1 state with ^ = 0 

is occupied by a fermion. Then, the Wfermion ground-state wavefunction is given by 

...,xjv) = det . (21) 

1<2 J<iV 

At hrst sight, one might think that the Slater determinants Eqs. flT^ and fl?Il) 

are not related to each other because {v^))“(a:)} and look completely different. 

Nevertheless, we will show that the following relation 

$^^‘'(a;i,...,a;^) = C^<l>"“(a;i,...,a;Ar) (22) 

holds for any N (even), where Cat is a constant independent of xi, ...,xn- This means 
that the many-body ground states of the uniform and SSD systems are identical to each 
other. The above relation can be proved by noting that the two Slater determinants 
become equivalent if there exists a linear transformation between {<p““(x)}()(rg^ and 
that is non-singular, i.e., the determinant of its transformation matrix is 
nonzero. For this purpose, however, the eigenfunctions in terms of the Gegenbauer 
polynomials ffT4|) are not so useful. In the following, we will consider a more direct 
treatment of ip^~^{x). Let us rewrite the intertwiner as 

A\ = — sin“^ (x) sin^ (x), (23) 

^ ax 

where we have introduced the variable x := ^x for simplicity. This leads us to the 
Rodrigues formula for the eigenfunctions Eq. (IT^ . 

iIj^~^{x) = sin”^’*'^ X \ -— ] sin^^“^ x. (24) 

^ \?>m.x ax 


This implies that the single-particle states can be described in terms of only 
trigonometric functions of x. Note that the explicit form depends on whether N and (3 
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are even or odd. Here, recall that N should be even for the PBC. We can thus write 
TV = 26* with a positive integer 6. For (3 = 2v = 0,1,6*), we get from Eq. fl2^ the 
following expansion: 

u^{.x) := (-1)^-'^ ^ ^ = 5^a2i.,2nsin2”(x), (25) 


where 


e-i 


0-2u,2n — 


k=n 


(2fc + l)(fc + 1) - I2{2l2 - 1) 

2{k-e){k + e) 


(26) 


with a 2 u, 2 e '■= 1- For (3 = 2u — 1 [u = 1,2,...,6), we also obtain another linearly- 
independent solution of odd parity as 


vp{x) := (-1) 


/TTX {2e + 2u-2)\ ,2(e-u)+l, 


= cos 


(40-1)! 


(x) a2u-i,2n-i sin^^ ^(x). 


-r2rr\x) 


(27) 


where 


e-i 

(^2u-l,2n-l = 

k=n 


k{ 2 k + l) - { 2 u-l){u-l) 

2{k-e){k + e) ^ 


(28) 


with 021 /-!, 20-1 := 1- 

Since i^) ^^^1 U 2 u{x) {v 2 u-iix)) are the same up to an 

overall constant, the only thing remaining to prove is to show that there is a linear 
transformation between the spaces of {m^(x), n^(x)} and {7^n“(^)}^=o 
invertible. From Eqs. fl25ll and fl27ll . one hnds the following relations: 

/ Uo{x)\ 

U2{x) 

\un{x)J 

( Vi{x) ^ 

■^ 3 ( 3 ;) 

\vn-i{x)J 

where matrix elements which are zero are left empty. These upper triangular matrices 
are non-singular and invertible since the diagonal element is nonzero for all (3. 
Thus, the linear space of the SSD eigenfunctions is equivalent to that of {sin^”(x)} 
and {cos(5:) sin^”“^(a:)}. According to the standard formula of the trigonometric 
functions [23], moreover, the linear spaces spanned by {sin^”(5:)} and {cos(x) sin^"“^(x)} 
are respectively equivalent to those by {cos(2nx)} and {sin(2nx)}. Hence, the linear 
space of {</?®®'^(x)}()ho can be spanned by the basis elements of {cos(2n5:)}^Zo ^^cl those 


/ 0-0,0 ^0,2 ■ 

0 - 2,2 • 

V 

/ ^1,1 ^1,3 • 

03,3 • 

V 


Oo,N \ 


( 1 ^ 

• 2 /~\ 

02 ,N 


Sin (x) 

on,n/ 


^sin'^(x) J 


(29) 


Ol,N-l \ ( cos(x) sin(x) \ 
Os^AT-i cos(5:) sin^(5;) 

on-i,n-i/ \cos(x) sin^“^(x)/ 


, (30) 
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of {sm(2nx)}^^i, indicating that the single-particle basis of the SSD problem is identical 
to that of the plane waves. 

Here, we should recall that the ^ = 0 state in the (3 = 0 sector is assumed to be 
empty of fermion, which implies that uq{x){o: ijjQ [x]) is not included in 
However, we can easily check Uq{x) = ®o.2n ^ cos(iVx) [21]. Thus, the 

removal of uq{x) from the SSD basis results in the linear space which is 

identical to that of the plane waves below the Fermi level in the uniform system. Thus, 
we have obtained the desired relation Eq. fl22l) 

If the total fermion number N is odd, we should take eigenfunctions in the row 
labeled by g = iV, which are antiperiodic in x. The proof proceeds along the same lines 
as that of the N =even case, and will not be given here. 

5. Single-particle states with positive energy 

The equivalence between the uniform ring and the SSD problem is established only for 
the ground states of many fermions. Here, we comment on single-particle excited states 
above the Fermi level, i.e., the positive energy > 0) solutions of the SSD problem. 
From the correspondence between the SSD and inverse sine-square problems, one can see 
that they translate into the eigenfunctions of the attractive 1 / sin^ potential problem. 
As noted in Appendix A, the solution of the 1 / sin^ potential problem with the attractive 
coupling is classified into two cases: (i) 0 < .^ < |(f)^ and (ii) < e The case 

(i) corresponds to Eq. (|H]) with 0 < < 1. For a non-integer value of (3, however, the 

solution of Eq. ([8]) is given by an analytic continuation of the Gegenbauer polynomial, 
which exhibits a singular behavior at the boundaries. Thus, the physical solution is 
permitted only in the limit of /5 —)■ 0 or 1, which is reduced to the plane-wave solution 
of /3 = 0 or 1. Recalling that the boundary condition for /3 = 0 is Neumann and that 
for = 1 is the Dirichlet, one may think of the solution for = 0 as the zero-energy 
limit of a scattering state, while that for (3 = 1 as the zero-energy limit of a bound 
state. For the case (ii), the attractive potential problem is ill-defined, where the energy 
of the single-particle ground state (in the sense of the potential problem) collapses to 
—oo [251 eq). Thus, the SSD problem does not have a proper single-particle excited 
state with positive energy, unlike the SSD for the lattice systems. We conclude that the 
SSD problem in continuous space is well defined only below the Fermi level. 

6. Geometric factor and lattice SSD 

In this section, we discuss the relation between the above continuous-space SSD and 
its lattice version, in the former of which a geometric factor associated with the inner 
product of the SSD eigenfunctions is of particular importance. So far, we have not 
defined the inner product of the SSD eigenfunctions. This is because we should pay a 
particular attention to the geometric factor in the continuous SSD, unlike the SSD for 
lattice systems. For instance, the standard definition of the inner product for the plane 
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waves fails to satisfy the orthogonality relation among (/9 = 1, 2,iV), though 

the linear space of them is identical to that of the plane waves. 

For the continuous SSD, the self-adjoint form of the differential equation tells us 
that an appropriate inner product for is 


((V', 


V-/3 




)) : = 


dx 


sm 


rX 


lN-0. 

I 


x) ij 


V' 


X 


= I 


N,0 


where the coefficient Jtv,^ is given by 


(31) 




{N + P-l)\ 


(2f)-l) (IV - f))! p - 1)!]: 


(32) 


the derivation of which is given in Appendix B The factor sin^ in the integrand in 
Eq. fl3l]) can be regarded as a geometric factor attributed to the spatial deformation. 
With this dehnition, the orthogonality of the single-particle wavefunctions for > 1 is 
correctly recovered. Here, we should mention that 'ipQ {x) cx cos{^x) for /3 = 0 is not 
normalizable with Eq. (ITIH . We will discuss this problem later. 

In order to see the role of the geometric factor, it is useful to consider the lattice 
version of the SSD problem, the Hamiltonian of which is defined as 




sm 


n=l 


T" I (4^+1 + 4+ic« 




sm 


71=1 


TT 

— In- 


4cn, (33) 


where cf and c„ are the creation and the annihilation operators, respectively, of a 
fermion at site n. The lattice spacing is set to unity so that the length of the chain 
L is the same as the total number of sites. We numerically solve the single-particle 
problem of T^iat and discuss its relation to the continuous SSD. Let Xmin) be the mth 
single-particle eigenfunction of "Hiat- The level index m runs from 0 to L — 1. The 
inner product between Xmip) and Xm'ip) is well-dehned and the following orthogonality 
holds: (Xm, Xm') ■= En=i = ^m, 7 n'• This suggests that, unlike Eq. ([M]), 

the lattice SSD eigenfunctions do not involve any geometric factor. On the other hand, it 
is also expected that a naive continuum limit of the lattice SSD eigenfunction Xm{n) may 
approach the eigenfunction of the continuous SSD problem 'ip^~^{x). These seemingly 
contradictory facts can be reconciled by comparing the inner product for the continuous 
SSD with that for the lattice ones. This leads to the following correspondence: 




N-j3 

/3 


(x) 



(34) 


with X ~ n, where sin(^n) can be viewed as a geometric correction to Xm- For 
convenience, we fix the relation between m and (3 as m = N — (3 so that the energy 
eigenvalues e(m) of the lattice eigenfunctions are arranged in ascending order. 

To verify Eq. (|3Tj) . we numerically diagonalize the lattice Hamiltonian (1331) for 
L = 101 with yU = 0 and obtain the single-particle eigenfunctions Xm{n). In this 
situation, we have 51 energy levels below the zero energy: m = 0, • • • , 50, where the 
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Figure 1. Eigenvalue spectrum of the lattice SSD model with L = 101 and 
/i = 0. The energy eigenvalues e(m) are arranged in ascending order. The Fermi 
level corresponds to e(50), which is exactly at zero energy. The solid curve indicates 
e(m) = —8.2 x (m — 49) (m — 50) which is inferred from the spectrum of the continuous 
SSD analytically obtained in section 3. 


exact zero energy state corresponds to m = 50. This can be seen in Fig. [H where the 
eigenvalue spectrum is shown as a function of m. Note that, because of the finite-L 
effect, the energy of the m = 49 state slightly deviates from the zero energy to the 
low-energy side. Thus, the many-body ground state of the lattice SSD problem can be 
obtained by Filing the levels m = 0, • • • ,49 with iV = 50 fermions. 

Before proceeding to the direct comparison of wavefunctions, we discuss the 
eigenvalue spectrum. Figure [T] shows the eigenvalue spectrum e(m) for the lattice SSD 
problem with L = 101 and /r = 0 in increasing order. In the hgure, we also plot the 
spectrum of the continuous SSD in the negative energy region: e(m) = —k(3{I3 — 1) with 
P = 50 — m. The overall coefficient k = 8.2 x 10“^ is determined so that the lattice 
and continuous results are on top of each other. Note that this value of k is clearly 
consistent with that of the continuous system k = ~ 9.67 x 10“^ for L = 101. As 

depicted in the hgure, we can see that the spectrum of the lattice SSD is well described 
by that of the continuous SSD in the entire energy range below the Fermi level. This 
should be contrasted to the fact that the linearization of the single-particle dispersion 
of the uniform lattice system, which is essential for taking the continuum limit, is valid 
only in the vicinity of the Fermi level. In this sense, the SSD for the lattice system 
efhciently smears out the lattice effect. 

Using the above correspondence between the lattice and continuous SSD, we can 
further explain the various characteristics of the spectrum around the Fermi energy 
(e = 0) reported in Ref. [8]. For instance, Eq. ([7]) for a hxed P{= N — m) leads to 

(35) 

with C := vr^/5(/5 — 1), which is in good agreement with the size dependence found in 
Ref. [8]. Note that this scaling relation contains no higher order correction of 0{L~^). 
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Figure 2. Comparison of single-particle wavefunctions between the lattice and 
continuous SSD systems. The chain length and the chemical potential are set to 
L = 101 and ^ = 0, respectively. In the corresponding continuous SSD, the total 
number of fermions is fixed as IV = 50. The lattice wavefunctions with the geometric 
correction are plotted as open circles for to = 0,1, 2, 25, while the continuum 

wavefunctions (' 07 v-m(®)) shown as solid lines. Note that the scaled position 
i := ^(n — i) is used for for comparison. 


Also from Eq. ([7]), we obtain the density of state (DOS) around the Fermi level as 


D(e) = 


djS 

de 


(lV 1 L I 
U; 2/1-1 “ 27r^ 


(36) 


for |e| 1/L^. For a sufficiently long chain, the DOS diverges with which 

is also consistent with the numerical result in Ref. [8]. For |e| l/D^, however, 

D ~ (^)^(=const), which is a cutoff bound in the vicinity of the Fermi level. 

Now, we turn to the comparison of the wavefunctions between the lattice and 
continuous SSD problems for L = 101 and N = 50. On the basis of Eq. flM)) . we define 
the lattice wavefunction with the geometric correction as 


:= sm(i) Xm(n) 


(37) 


with X = — |). The spatial prohles of for m = 0,1, 2, 25 are shown in Fig. [2l 

The corresponding wavefunctions of the continuous SSD, '0]^_^(a;), for m = 0,1,2,25 
are also plotted as solid lines, where the overall normalization is adjusted so that 
and almost on top of each other. One can clearly see that the single-particle 

ground-state {m = N — (3 = 0) is localized at the center of the chain and ■^^'^(n) tends 
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to delocalize as m increases. It is remarkable that reproduces 

even when m = 25, which exhibits a rapid oscillation as a function of x. Here, we note 
that the geometric correction becomes effective near the boundaries of the system. If 
we omit this correction, the lattice wavefunction signihcantly deviates from 
as X ^ 0 or TT. We thus conclude that the correspondence fIMl) between lattice and 
continuous SSD is well established not only at the level of energy eigenvalues, but also 
at the level of individual eigenfunctions. This implies that the continuum limit of the 
lattice is taken efficiently in the SSD system, which is consistent with the fact that the 
leading hnite-size dependence of the excitation energy is proportional to 1/T^, instead 
of 1/L which is universal in CFT. 



Figure 3. The lattice SSD wavefunction at the Fermi level and geometric correction 
for the case of L = 101 and /i = 0. The blue broken line with open circles indicates 
X 5 o(n) as a function of the scaled position x = ^{n — The red solid line with open 
circles represents = sin(i)x 5 o(n). 

The geometric correction signihcantly affects the wavefunctions near the Fermi 
level, which is extended over the entire region of the chain. In Fig. [3l we show the 
lattice wavefunction X 5 o(^) for L = 101 with /i = 0, where m = 50 corresponds to the 
zero energy state just on the Fermi level. In the hgure, we also plot oc '0o°(^) 

for comparison. Then, it is observed that the amplitude of the lattice wavefunction 
X 5 o(n) increases toward the end sites n = 1 and L, while ' 05 o^(^) ~ 'i/^o°(^) retains a 
plane-wave behavior. As was mentioned before, 'ipQ {x) is not normalizable with respect 
to the inner product Eq. fl3lli since '0^ {x)/sm{Nx) ~ cot(A^x) diverges at a; = 0 and tt. 
This divergence comes from the fact that i/jq {x) has a hnite amplitude at the edges due 
to the Neumann boundary conditions. The result of Fig. [3] implies that the divergence 
of the inner product flTIli for , which originates solely from the geometric factor, can 
be regularized by a lattice cutoff. 
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7. summary and discussions 

In this paper, we have studied the sine-square deformation (SSD) of the free-fermion 
problem in one-dimensional continuous space. We formulated the single-particle problem 
of the SSD system as the inverse of the 1/sin^ potential problem, where supersymmetric 
(SUSY) quantum mechanics is very effective. Exploiting both SUSY and the shape 
invariance, we explicitly obtained the complete set of the single-particle eigenfunctions 
of the SSD problem, and showed that its linear space is identical to that of the uniform 
periodic (ant ip er iodic) ring. Accordingly, we have shown that the many-body ground 
state of the system with SSD is the same as that of the uniform system with periodic 
(antiperiodic) boundary conditions. It would be interesting to generalize our theory to 
other class of superpotentials with the shape invariance. In addition, we have found that 
various properties of the lattice SSD problem are well explained by the analytic results 
for the SSD in the continuous space. In particular, the eigenfunctions of the lattice SSD 
problem can be reproduced from those of the continuous SSD by taking into account 
the geometric correction. 

In the previous studies of the lattice systems that are reducible to free fermions [11 
El [12], the ground-state correspondence was proved using the plane-wave basis. Also, the 
correspondence was explained in the context of GET [5|, where the SSD Hamiltonian 
is expressed as Lq — (L+i -|- L_i)/2 and both the Virasoro generators L± annihilate 
the vacuum state |0). The present study clearly provides a complementary approach 
to the SSD problem, and various aspects of the SSD eigenfunctions can be clarified 
with the help of SUSY quantum mechanics. For example, we have demonstrated in 
the SSD system that the 1/L^ correction fl35l) of the excitation energy emerges with no 
higher order terms, instead of the universal 1/L-size dependence in GET. In this paper, 
nevertheless, SUSY quantum mechanics was just used to systematically construct the 
single-particle eigenfunctions of the SSD problem, and thus the role of the SUSY may not 
be clear at the level of held theory. Moreover, a very recent paper m pointed out that an 
inhnite circumference limit of GET can be obtained by adopting Lq —(T_|_i-|-L_i)/2 as the 
Hamiltonian instead of Lq. Such singular behavior may be related to the normalizability 
of the zero energy state satisfying the Neumann boundary conditions that appear in 
the present continuous SSD. Therefore, it would be interesting to explore the relation 
between SUSY and GET behind the SSD problems. 
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Appendix A. Formal solutions for positive energy 

Appendix A.l. 0 < .^ < ^ 

For Eq. (|8]), ^ becomes positive in 0 < /3 < 1, which could be a candidate of a scattering 
state with positive energy. Taking account of a non-integer valued fd, we introduce 
another parameterization jd = a + 1/2, which yields 

! - - (B’ (•■ - y ■ 

We then obtain the linearly-independent solutions of Eq. ([8]) for a hxed pi = (f as 
u{x) = (sinx)^/2+“cosxF Q + f ^ cos^x^ (A.2) 

n(x) = (sinx)^/2+“F Q + | ^ f cos^x^ , (A.3) 

where F is Gauss’s hypergeometric function. Note that, if a is a positive half integer, 
Eqs. flA.2p and flA.Sp reduce to the Gegenbauer polynomials. Otherwise, these are not 
a polynomials in cos^ x. 

Since Eqs. flA.2p and (lA.3p with \a\ < 1/2 are not divergent in 0 < x < L, one 
might expect that the scattering state could be described by a linear combination of 
them. However, the derivative of u{x) and v{x) with respect to x diverges at x = 0 
and 1 when |a| < 1/2, as far as q is fixed at an integer [27]. Thus, there is no physical 
solution when 0 < e < lijY- Of course, in the limit of a —)■ ±1/2, u{x) and v{x) 
approach to the trigonometric functions that are solutions for fd = 0 and fd = 1. 


Appendix A.2. ^ 


In the context of the 1/sin^ potential problem, the coupling corresponding to ^ | (f 

is so strongly attractive that the single-particle ground state of is unstable, i.e., p, 
is not bounded from below. Formal solutions of Eq. dHj) for a certain hxed chemical 
potential p are also given by Gauss’s hypergeometric function 


u{x) 


v{x) 


(sinx)^/2+*“cosxF Q + 
(sinx)^/2+“FQ + *|- 


.a pi .a pi 2 ~ 

i -, —V I —I—; cos x 

2 2 4 2 2 2 


p 1 a p 1 2 

—, T ± *——; cos X 
2 4 2 2 2 


(A.4) 

(A.5) 


where we have adopted the parameterization ^ (o^ + i)- However, Eqs. flA.4p and 

flA.5ll exhibit singular behavior in the vicinities of x = 0 and 1, rehecting the unstable 
ground state. 


Appendix B. Derivation of in Eq. (1311) 

In this Appendix, we present a detailed derivation of the coefhcient appearing in 
Eq. fl3T]l . Gonsider the squared norm of '^’^{x), where fd is assumed to be an integer and 
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/9 > 1. Changing the variable from x to z := cos ^x, we have 

{{'^ 0 , i’D) = [ dz[C^{z)]^l-zy-^. 


TT 


(B.l) 


'-1 


To evaluate the above integral, we expand the Gegenbauer polynomials in terms of 
{n = 0,1, The results depend on whether I is even or odd: 




fc =0 


/3-1 


fc =0 


/3-1 


which can be proved using the formula obtained by Gegenbauer himself 
using the orthogonality of C^~^{z)\ 


/: 

we arrive at 


dzC, [z)C, (z)(l z) - + V, 


2^ 2/5 (A7j)2 ^ 


(B.2) 

Then, 

(B.3) 

(B.4) 


The coefficient can be read off from the above formula by replacing £ with N — (3. 
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